Genome-resolved metagenomics of Venice Lagoon surface sediment bacteria reveals high biosynthetic potential and metabolic plasticity as successful strategies in an impacted environment.

Bacteria living in sediments play essential roles in marine ecosystems and deeper insights into the ecology and biogeochemistry of these largely unexplored organisms can be obtained from 'omics' approaches. Here, we characterized metagenome-assembled-genomes (MAGs) from the surface sediment microbes of the Venice Lagoon (northern Adriatic Sea) in distinct sub-basins exposed to various natural and anthropogenic pressures. MAGs were explored for biodiversity, major marine metabolic processes, anthropogenic activity-related functions, adaptations at the microscale, and biosynthetic gene clusters. Starting from 126 MAGs, a non-redundant dataset of 58 was compiled, the majority of which (35) belonged to (Alpha- and Gamma-) Proteobacteria. Within the broad microbial metabolic repertoire (including C, N, and S metabolisms) the potential to live without oxygen emerged as one of the most important features. Mixotrophy was also found as a successful lifestyle. Cluster analysis showed that different MAGs encoded the same metabolic patterns (e.g., C fixation, sulfate oxidation) thus suggesting metabolic redundancy. Antibiotic and toxic compounds resistance genes were coupled, a condition that could promote the spreading of these genetic traits. MAGs showed a high biosynthetic potential related to antimicrobial and biotechnological classes and to organism defense and interactions as well as adaptive strategies for micronutrient uptake and cellular detoxification. Our results highlighted that bacteria living in an impacted environment, such as the surface sediments of the Venice Lagoon, may benefit from metabolic plasticity as well as from the synthesis of a wide array of secondary metabolites, promoting ecosystem resilience and stability toward environmental pressures. Supplementary Information: The online version contains supplementary material available at 10.1007/s42995-023-00192-z.

Annual recurrence of prokaryotic climax communities in shallow waters of the North Mediterranean.

In temperate coastal environments, wide fluctuations of biotic and abiotic factors drive microbiome dynamics. To link recurrent ecological patterns with planktonic microbial communities, we analysed a monthly-sampled 3-year time series of 16S rRNA amplicon sequencing data, alongside environmental variables, collected at two stations in the northern Adriatic Sea. Time series multivariate analyses allowed us to identify three stable, mature communities (climaxes), whose recurrence was mainly driven by changes in photoperiod and temperature. Mixotrophs (e.g., Ca. Nitrosopumilus, SUP05 clade, and Marine Group II) thrived under oligotrophic, low-light conditions, whereas copiotrophs (e.g., NS4 and NS5 clades) bloomed at higher temperatures and substrate availability. The early spring climax was characterised by a more diverse set of amplicon sequence variants, including copiotrophs associated with phytoplankton-derived organic matter degradation, and photo-auto/heterotrophic organisms (e.g., Synechococcus sp., Roseobacter clade), whose rhythmicity was linked to photoperiod lengthening. Through the identification of recurrent climax assemblages, we begin to delineate a typology of ecosystem based on microbiome composition and functionality, allowing for the intercomparison of microbial assemblages among different biomes, a still underachieved goal in the omics era.

The Impact of MOSE (Experimental Electromechanical Module) Flood Barriers on Microphytobenthic Community of the Venice Lagoon.

MOSE is a system of mobile gates engineered to temporarily isolate the Venice Lagoon from the Adriatic Sea and to protect the city from flooding during extreme high tides. Within the framework of the Venezia2021 program, we conducted two enclosure experiments in July 2019 (over 48 h) and October 2020 (over 28 h) by means of 18 mesocosms, in order to simulate the structural alterations that microphytobenthos (MPB) assemblages might encounter when the MOSE system is operational. The reduced hydrodynamics inside the mesocosms favored the deposition of organic matter and the sinking of cells from the water column towards the sediment. Consequently, MPB abundances increased over the course of both experiments and significant changes in the taxonomic composition of the community were recorded. Species richness increased in summer while it slightly decreased in autumn, this latter due to the increase in relative abundances of taxa favored by high organic loads and fine grain size. By coupling classical taxonomy with 18S rRNA gene metabarcoding we were able to obtain a comprehensive view of the whole community potential, highlighting the complementarity of these two approaches in ecological studies. Changes in the structure of MPB could affect sediment biostabilization, water turbidity and lagoon primary production.

Improving environmental monitoring of Vibrionaceae in coastal ecosystems through 16S rRNA gene amplicon sequencing.

The Vibrionaceae family groups genetically and metabolically diverse bacteria thriving in all marine environments. Despite often representing a minor fraction of bacterial assemblages, members of this family can exploit a wide variety of nutritional sources, which makes them important players in biogeochemical dynamics. Furthermore, several Vibrionaceae species are well-known pathogens, posing a threat to human and animal health. Here, we applied the phylogenetic placement coupled with a consensus-based approach using 16S rRNA gene amplicon sequencing, aiming to reach a reliable and fine-level Vibrionaceae characterization and identify the dynamics of blooming, ecologically important, and potentially pathogenic species in different sites of the northern Adriatic Sea. Water samples were collected monthly at a Long-Term Ecological Research network site from 2018 to 2021, and in spring and summer of 2019 and 2020 at two sites affected by depurated sewage discharge. The 41 identified Vibrionaceae species represented generally below 1% of the sampled communities; blooms (up to ~ 11%) mainly formed by Vibrio chagasii and Vibrio owensii occurred in summer, linked to increasing temperature and particulate matter concentration. Pathogenic species such as Vibrio anguilllarum, Vibrio tapetis, and Photobacterium damselae were found in low abundance. Depuration plant samples were characterized by a lower abundance and diversity of Vibrionaceae species compared to seawater, highlighting that Vibrionaceae dynamics at sea are unlikely to be related to wastewater inputs. Our work represents a further step to improve the molecular approach based on short reads, toward a shared, updated, and curated phylogeny of the Vibrionaceae family.

Linking lifestyle and foraging strategies of marine bacteria: selfish behaviour of particle-attached bacteria in the northern Adriatic Sea.

Microbe-mediated enzymatic hydrolysis of organic matter entails the production of hydrolysate, the recovery of which may be more or less efficient. The selfish uptake mechanism, recently discovered, allows microbes to hydrolyze polysaccharides and take up large oligomers, which are then degraded in the periplasmic space. By minimizing the hydrolysate loss, selfish behaviour may be profitable for free-living cells dwelling in a patchy substrate landscape. However, selfish uptake seems to be tailored to algal-derived polysaccharides, abundant in organic particles, suggesting that particle-attached microbes may use this strategy. We tracked selfish polysaccharides uptake in surface microbial communities of the northeastern Mediterranean Sea, linking the occurrence of this processing mode with microbial lifestyle. Additionally, we set up fluorescently labelled polysaccharides incubations supplying phytodetritus to investigate a 'pioneer' scenario for particle-attached microbes. Under both conditions, selfish behaviour was almost exclusively carried out by particle-attached microbes, suggesting that this mechanism may represent an advantage in the race for particle exploitation. Our findings shed light on the selfish potential of particle-attached microbes, suggesting multifaceted foraging strategies exerted by particle colonizers.

High-Throughput Sequencing of Environmental DNA as a Tool for Monitoring Eukaryotic Communities and Potential Pathogens in a Coastal Upwelling Ecosystem.

The marine environment includes diverse microeukaryotic organisms that play important functional roles in the ecosystem. With molecular approaches, eukaryotic taxonomy has been improved, complementing classical analysis. In this study, DNA metabarcoding was performed to describe putative pathogenic eukaryotic microorganisms in sediment and marine water fractions collected in Galicia (NW Spain) from 2016 to 2018. The composition of eukaryotic communities was distinct between sediment and water fractions. Protists were the most diverse group, with the clade TSAR (Stramenopiles, Alveolata, Rhizaria, and Telonemida) as the primary representative organisms in the environment. Harmful algae and invasive species were frequently detected. Potential pathogens, invasive pathogenic organisms as well as the causative agents of harmful phytoplanktonic blooms were identified in this marine ecosystem. Most of the identified pathogens have a crucial impact on the aquacultural sector or affect to relevant species in the marine ecosystem, such as diatoms. Moreover, pathogens with medical and veterinary importance worldwide were also found, as well as pathogens that affect diatoms. The evaluation of the health of a marine ecosystem that directly affects the aquacultural sector with a zoonotic concern was performed with the metabarcoding assay.

Environmental DNA assessment of airborne plant and fungal seasonal diversity.

Environmental DNA (eDNA) metabarcoding and metagenomics analyses can improve taxonomic resolution in biodiversity studies. Only recently, these techniques have been applied in aerobiology, to target bacteria, fungi and plants in airborne samples. Here, we present a nine-month aerobiological study applying eDNA metabarcoding in which we analyzed simultaneously airborne diversity and variation of fungi and plants across five locations in North and Central Italy. We correlated species composition with the ecological characteristics of the sites and the seasons. The most abundant taxa among all sites and seasons were the fungal genera Cladosporium, Alternaria, and Epicoccum and the plant genera Brassica, Corylus, Cupressus and Linum, the latter being much more variable among sites. PERMANOVA and indicator species analyses showed that the plant diversity from air samples is significantly correlated with seasons, while that of fungi varied according to the interaction between seasons and sites. The results consolidate the performance of a new eDNA metabarcoding pipeline for the simultaneous amplification and analysis of airborne plant and fungal particles. They also highlight the promising complementarity of this approach with more traditional biomonitoring frameworks and routine reports of air quality provided by environmental agencies.

Prokaryotic Response to Phytodetritus-Derived Organic Material in Epi- and Mesopelagic Antarctic Waters.

Particulate organic matter (POM) export represents the underlying principle of the biological carbon pump, driving the carbon flux from the sunlit to the dark ocean. The efficiency of this process is tightly linked to the prokaryotic community, as >70% of POM respiration is carried out by particle-associated prokaryotes. In the Ross Sea, one of the most productive areas of the Southern Ocean, up to 50% of the surface primary production is exported to the mesopelagic ocean as POM. Recent evidence suggests that a significant fraction of the POM in this area is composed of intact phytoplankton cells. During austral summer 2017, we set up bottle enrichment experiments in which we amended free-living surface and deep prokaryotic communities with organic matter pools generated from native microplankton, mimicking the particle export that may derive from mild (1 µg of Chlorophyll a L-1) and intense (10 µg of Chlorophyll a L-1) phytoplankton bloom. Over a course of 4 days, we followed free-living and particle-attached prokaryotes' abundance, the degradation rates of polysaccharides, proteins and lipids, heterotrophic production as well as inorganic carbon utilization and prokaryotic community structure dynamics. Our results showed that several rare or undetected taxa in the initial community became dominant during the time course of the incubations and that different phytodetritus-derived organic matter sources induced specific changes in microbial communities, selecting for peculiar degradation and utilization processes spectra. Moreover, the features of the supplied detritus (in terms of microplankton taxa composition) determined different colonization dynamics and organic matter processing modes. Our study provides insights into the mechanisms underlying the prokaryotic utilization of phytodetritus, a significant pool of organic matter in the dark ocean.

PLANiTS: a curated sequence reference dataset for plant ITS DNA metabarcoding.

DNA metabarcoding combines DNA barcoding with high-throughput sequencing to identify different taxa within environmental communities. The ITS has already been proposed and widely used as universal barcode marker for plants, but a comprehensive, updated and accurate reference dataset of plant ITS sequences has not been available so far. Here, we constructed reference datasets of Viridiplantae ITS1, ITS2 and entire ITS sequences including both Chlorophyta and Streptophyta. The sequences were retrieved from NCBI, and the ITS region was extracted. The sequences underwent identity check to remove misidentified records and were clustered at 99% identity to reduce redundancy and computational effort. For this step, we developed a script called 'better clustering for QIIME' (bc4q) to ensure that the representative sequences are chosen according to the composition of the cluster at a different taxonomic level. The three datasets obtained with the bc4q script are PLANiTS1 (100 224 sequences), PLANiTS2 (96 771 sequences) and PLANiTS (97 550 sequences), and all are pre-formatted for QIIME, being this the most used bioinformatic pipeline for metabarcoding analysis. Being curated and updated reference databases, PLANiTS1, PLANiTS2 and PLANiTS are proposed as a reliable, pivotal first step for a general standardization of plant DNA metabarcoding studies. The bc4q script is presented as a new tool useful in each research dealing with sequences clustering. Database URL: https://github.com/apallavicini/bc4q; https://github.com/apallavicini/PLANiTS.

Graphene-based materials do not impair physiology, gene expression and growth dynamics of the aeroterrestrial microalga Trebouxia gelatinosa.

The effects of two graphene-based materials (GBMs), few-layers graphene (FLG) and graphene oxide (GO), were studied in the aeroterrestrial green microalga Trebouxia gelatinosa. Algae were subjected to short- and long-term exposure to GBMs at 0.01, 1 and 50 µg mL - 1. GBMs internalization after short-term exposures was investigated with confocal microscopy, Raman spectroscopy and TEM. Potential negative effects of GBMs, compared to the oxidative stress induced by H2O2, were verified by analyzing chlorophyl a fluorescence (ChlaF), expression of stress-related genes and membrane integrity. Effects of up to 4-week-long exposures were assessed analyzing growth dynamics, ChlaF and photosynthetic pigments. GBMs were not observed in cells but FLG was detected at the interface between the cell wall and plasma membrane, whereas GO was observed adherent to the external wall surface. FLG caused the down-regulation of the HSP70-1 gene, with the protein levels remaining stable, whereas GO had no effect. In comparison, H2O2 produced dose- and time-dependent effects on ChlaF, gene expression and HSP70 protein level. Long-term exposures to GBMs did not affect growth dynamics, ChlaF or photosynthetic pigment contents, indicating that the few observed short-term effects were not dangerous on the long-term. Results suggest that interactions between FLG and plasma membrane were harmless, activating a down-regulation of the HSP70-1 gene similar to that induced by H2O2. Our work shows that studying GBMs effects on non-model organisms is important since the results of model green microalgae are not representative of the whole taxonomic group.

The Rpv3-3 Haplotype and Stilbenoid Induction Mediate Downy Mildew Resistance in a Grapevine Interspecific Population.

The development of new resistant varieties to the oomycete Plasmopara viticola (Berk.& Curt) is a promising way to combat downy mildew (DM), one of the major diseases threatening the cultivated grapevine (Vitis vinifera L.). Taking advantage of a segregating population derived from "Merzling" (a mid-resistant hybrid) and "Teroldego" (a susceptible landrace), 136 F1 individuals were characterized by combining genetic, phenotypic, and gene expression data to elucidate the genetic basis of DM resistance and polyphenol biosynthesis upon P. viticola infection. An improved consensus linkage map was obtained by scoring 192 microsatellite markers. The progeny were screened for DM resistance and production of 42 polyphenols. QTL mapping showed that DM resistance is associated with the herein named Rpv3-3 specific haplotype and it identified 46 novel metabolic QTLs linked to 30 phenolics-related parameters. A list of the 95 most relevant candidate genes was generated by specifically exploring the stilbenoid-associated QTLs. Expression analysis of 11 genes in Rpv3-3 +/- genotypes displaying disparity in DM resistance level and stilbenoid accumulation revealed significant new candidates for the genetic control of stilbenoid biosynthesis and oligomerization. These overall findings emphasized that DM resistance is likely mediated by the major Rpv3-3 haplotype and stilbenoid induction.

Relation between water status and desiccation-affected genes in the lichen photobiont Trebouxia gelatinosa.

The relation between water status and expression profiles of desiccation -related genes has been studied in the desiccation tolerant (DT) aeroterrestrial green microalga Trebouxia gelatinosa, a common lichen photobiont. Algal colonies were desiccated in controlled conditions and during desiccation water content (WC) and water potential (Ψ) were measured to find the turgor loss point (Ψtlp). Quantitative real-time PCR was performed to measure the expression of ten genes related to photosynthesis, antioxidant defense, expansins, heat shock proteins (HSPs), and desiccation related proteins in algal colonies collected during desiccation when still at full turgor (WC > 6 g H2O g-1 dry weight), immediately before and after Ψtlp (-4 MPa; WC ∼ 1 g H2O g-1 dry weight) and before and after complete desiccation (WC < 0.01 g H2O g-1 dry weight), quantifying the HSP70 protein levels by immunodetection. Our analysis showed that the expression of eight out of ten genes changed immediately before and after Ψtlp. Interestingly, the expression of five out of ten genes changed also before complete desiccation, i.e. between 0.2 and 0.01 g H2O g-1 dry weight. However, the HSP70 protein levels were not affected by changes in water status. The study provides new evidences of the link between the loss of turgor and the expression of genes related to the desiccation tolerance of T. gelatinosa, suggesting the former as a signal triggering inducible mechanisms.

DNA metabarcoding uncovers fungal diversity of mixed airborne samples in Italy.

Fungal spores and mycelium fragments are particles which become and remain airborne and have been subjects of aerobiological studies. The presence and the abundance of taxa in aerobiological samples can be very variable and impaired by changeable climatic conditions. Because many fungi produce mycotoxins and both their mycelium fragments and spores are potential allergens, monitoring the presence of these taxa is of key importance. So far data on exposure and sensitization to fungal allergens are mainly based on the assessment of few, easily identifiable taxa and focused only on certain environments. The microscopic method used to analyze aerobiological samples and the inconspicuous fungal characters do not allow a in depth taxonomical identification. Here, we present a first assessment of fungal diversity from airborne samples using a DNA metabarcoding analysis. The nuclear ITS2 region was selected as barcode to catch fungal diversity in mixed airborne samples gathered during two weeks in four sites of North-Eastern and Central Italy. We assessed the taxonomic composition and diversity within and among the sampled sites and compared the molecular data with those obtained by traditional microscopy. The molecular analyses provide a tenfold more comprehensive determination of the taxa than the traditional morphological inspections. Our results prove that the metabarcoding analysis is a promising approach to increases quality and sensitivity of the aerobiological monitoring. The laboratory and bioinformatic workflow implemented here is now suitable for routine, high-throughput, regional analyses of airborne fungi.

ITS2 metabarcoding analysis complements lichen mycobiome diversity data.

Lichen thalli harbor complex fungal communities (mycobiomes) of species with divergent trophic and ecological strategies. The complexity and diversity of lichen mycobiomes are still largely unknown, despite surveys combining culture-based methods and high-throughput sequencing (HTS). The results of such surveys are strongly influenced by the barcode locus chosen, its sensitivity in discriminating taxa, and the depth to which public sequence repositories cover the phylogenetic spectrum of fungi. Here, we use HTS of the internal transcribed spacer 2 (ITS2) to assess the taxonomic composition and diversity of a well-characterized, alpine rock lichen community that includes thalli symptomatically infected by lichenicolous fungi as well as asymptomatic thalli. Taxa belonging to the order Chaetothyriales are the major components of the observed lichen mycobiomes. We predict sequences representative of lichenicolous fungi characterized morphologically and assess their asymptomatic presence in lichen thalli. We demonstrated the limitations of metabarcoding in fungi and show how the estimation of species diversity widely differs when ITS1 or ITS2 are used as barcode, and particularly biases the detection of Basidiomycota. The complementary analysis of both ITS1 and ITS2 loci is therefore required to reliably estimate the diversity of lichen mycobiomes.

Integrated QTL detection for key breeding traits in multiple peach progenies.

BACKGROUND: Peach (Prunus persica (L.) Batsch) is a major temperate fruit crop with an intense breeding activity. Breeding is facilitated by knowledge of the inheritance of the key traits that are often of a quantitative nature. QTLs have traditionally been studied using the phenotype of a single progeny (usually a full-sib progeny) and the correlation with a set of markers covering its genome. This approach has allowed the identification of various genes and QTLs but is limited by the small numbers of individuals used and by the narrow transect of the variability analyzed. In this article we propose the use of a multi-progeny mapping strategy that used pedigree information and Bayesian approaches that supports a more precise and complete survey of the available genetic variability. RESULTS: Seven key agronomic characters (data from 1 to 3 years) were analyzed in 18 progenies from crosses between occidental commercial genotypes and various exotic lines including accessions of other Prunus species. A total of 1467 plants from these progenies were genotyped with a 9 k SNP array. Forty-seven QTLs were identified, 22 coinciding with major genes and QTLs that have been consistently found in the same populations when studied individually and 25 were new. A substantial part of the QTLs observed (47%) would not have been detected in crosses between only commercial materials, showing the high value of exotic lines as a source of novel alleles for the commercial gene pool. Our strategy also provided estimations on the narrow sense heritability of each character, and the estimation of the QTL genotypes of each parent for the different QTLs and their breeding value. CONCLUSIONS: The integrated strategy used provides a broader and more accurate picture of the variability available for peach breeding with the identification of many new QTLs, information on the sources of the alleles of interest and the breeding values of the potential donors of such valuable alleles. These results are first-hand information for breeders and a step forward towards the implementation of DNA-informed strategies to facilitate selection of new cultivars with improved productivity and quality.

New features of desiccation tolerance in the lichen photobiont Trebouxia gelatinosa are revealed by a transcriptomic approach.

Trebouxia is the most common lichen-forming genus of aero-terrestrial green algae and all its species are desiccation tolerant (DT). The molecular bases of this remarkable adaptation are, however, still largely unknown. We applied a transcriptomic approach to a common member of the genus, T. gelatinosa, to investigate the alteration of gene expression occurring after dehydration and subsequent rehydration in comparison to cells kept constantly hydrated. We sequenced, de novo assembled and annotated the transcriptome of axenically cultured T. gelatinosa by using Illumina sequencing technology. We tracked the expression profiles of over 13,000 protein-coding transcripts. During the dehydration/rehydration cycle c. 92 % of the total protein-coding transcripts displayed a stable expression, suggesting that the desiccation tolerance of T. gelatinosa mostly relies on constitutive mechanisms. Dehydration and rehydration affected mainly the gene expression for components of the photosynthetic apparatus, the ROS-scavenging system, Heat Shock Proteins, aquaporins, expansins, and desiccation related proteins (DRPs), which are highly diversified in T. gelatinosa, whereas Late Embryogenesis Abundant Proteins were not affected. Only some of these phenomena were previously observed in other DT green algae, bryophytes and resurrection plants, other traits being distinctive of T. gelatinosa, and perhaps related to its symbiotic lifestyle. Finally, the phylogenetic inference extended to DRPs of other chlorophytes, embryophytes and bacteria clearly pointed out that DRPs of chlorophytes are not orthologous to those of embryophytes: some of them were likely acquired through horizontal gene transfer from extremophile bacteria which live in symbiosis within the lichen thallus.

Whole-Genome Analysis of Diversity and SNP-Major Gene Association in Peach Germplasm.

Peach was domesticated in China more than four millennia ago and from there it spread world-wide. Since the middle of the last century, peach breeding programs have been very dynamic generating hundreds of new commercial varieties, however, in most cases such varieties derive from a limited collection of parental lines (founders). This is one reason for the observed low levels of variability of the commercial gene pool, implying that knowledge of the extent and distribution of genetic variability in peach is critical to allow the choice of adequate parents to confer enhanced productivity, adaptation and quality to improved varieties. With this aim we genotyped 1,580 peach accessions (including a few closely related Prunus species) maintained and phenotyped in five germplasm collections (four European and one Chinese) with the International Peach SNP Consortium 9K SNP peach array. The study of population structure revealed the subdivision of the panel in three main populations, one mainly made up of Occidental varieties from breeding programs (POP1OCB), one of Occidental landraces (POP2OCT) and the third of Oriental accessions (POP3OR). Analysis of linkage disequilibrium (LD) identified differential patterns of genome-wide LD blocks in each of the populations. Phenotypic data for seven monogenic traits were integrated in a genome-wide association study (GWAS). The significantly associated SNPs were always in the regions predicted by linkage analysis, forming haplotypes of markers. These diagnostic haplotypes could be used for marker-assisted selection (MAS) in modern breeding programs.

Accuracy and responses of genomic selection on key traits in apple breeding.

The application of genomic selection in fruit tree crops is expected to enhance breeding efficiency by increasing prediction accuracy, increasing selection intensity and decreasing generation interval. The objectives of this study were to assess the accuracy of prediction and selection response in commercial apple breeding programmes for key traits. The training population comprised 977 individuals derived from 20 pedigreed full-sib families. Historic phenotypic data were available on 10 traits related to productivity and fruit external appearance and genotypic data for 7829 SNPs obtained with an Illumina 20K SNP array. From these data, a genome-wide prediction model was built and subsequently used to calculate genomic breeding values of five application full-sib families. The application families had genotypes at 364 SNPs from a dedicated 512 SNP array, and these genotypic data were extended to the high-density level by imputation. These five families were phenotyped for 1 year and their phenotypes were compared to the predicted breeding values. Accuracy of genomic prediction across the 10 traits reached a maximum value of 0.5 and had a median value of 0.19. The accuracies were strongly affected by the phenotypic distribution and heritability of traits. In the largest family, significant selection response was observed for traits with high heritability and symmetric phenotypic distribution. Traits that showed non-significant response often had reduced and skewed phenotypic variation or low heritability. Among the five application families the accuracies were uncorrelated to the degree of relatedness to the training population. The results underline the potential of genomic prediction to accelerate breeding progress in outbred fruit tree crops that still need to overcome long generation intervals and extensive phenotyping costs.

Development and validation of a 20K single nucleotide polymorphism (SNP) whole genome genotyping array for apple (Malus × domestica Borkh).

High-density SNP arrays for genome-wide assessment of allelic variation have made high resolution genetic characterization of crop germplasm feasible. A medium density array for apple, the IRSC 8K SNP array, has been successfully developed and used for screens of bi-parental populations. However, the number of robust and well-distributed markers contained on this array was not sufficient to perform genome-wide association analyses in wider germplasm sets, or Pedigree-Based Analysis at high precision, because of rapid decay of linkage disequilibrium. We describe the development of an Illumina Infinium array targeting 20K SNPs. The SNPs were predicted from re-sequencing data derived from the genomes of 13 Malus × domestica apple cultivars and one accession belonging to a crab apple species (M. micromalus). A pipeline for SNP selection was devised that avoided the pitfalls associated with the inclusion of paralogous sequence variants, supported the construction of robust multi-allelic SNP haploblocks and selected up to 11 entries within narrow genomic regions of ±5 kb, termed focal points (FPs). Broad genome coverage was attained by placing FPs at 1 cM intervals on a consensus genetic map, complementing them with FPs to enrich the ends of each of the chromosomes, and by bridging physical intervals greater than 400 Kbps. The selection also included ∼3.7K validated SNPs from the IRSC 8K array. The array has already been used in other studies where ∼15.8K SNP markers were mapped with an average of ∼6.8K SNPs per full-sib family. The newly developed array with its high density of polymorphic validated SNPs is expected to be of great utility for Pedigree-Based Analysis and Genomic Selection. It will also be a valuable tool to help dissect the genetic mechanisms controlling important fruit quality traits, and to aid the identification of marker-trait associations suitable for the application of Marker Assisted Selection in apple breeding programs.

Evaluation of SNP Data from the Malus Infinium Array Identifies Challenges for Genetic Analysis of Complex Genomes of Polyploid Origin.

High throughput arrays for the simultaneous genotyping of thousands of single-nucleotide polymorphisms (SNPs) have made the rapid genetic characterisation of plant genomes and the development of saturated linkage maps a realistic prospect for many plant species of agronomic importance. However, the correct calling of SNP genotypes in divergent polyploid genomes using array technology can be problematic due to paralogy, and to divergence in probe sequences causing changes in probe binding efficiencies. An Illumina Infinium II whole-genome genotyping array was recently developed for the cultivated apple and used to develop a molecular linkage map for an apple rootstock progeny (M432), but a large proportion of segregating SNPs were not mapped in the progeny, due to unexpected genotype clustering patterns. To investigate the causes of this unexpected clustering we performed BLAST analysis of all probe sequences against the 'Golden Delicious' genome sequence and discovered evidence for paralogous annealing sites and probe sequence divergence for a high proportion of probes contained on the array. Following visual re-evaluation of the genotyping data generated for 8,788 SNPs for the M432 progeny using the array, we manually re-scored genotypes at 818 loci and mapped a further 797 markers to the M432 linkage map. The newly mapped markers included the majority of those that could not be mapped previously, as well as loci that were previously scored as monomorphic, but which segregated due to divergence leading to heterozygosity in probe annealing sites. An evaluation of the 8,788 probes in a diverse collection of Malus germplasm showed that more than half the probes returned genotype clustering patterns that were difficult or impossible to interpret reliably, highlighting implications for the use of the array in genome-wide association studies.